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Abstract. We consider a small ensemble of Bose atoms in a ring optical lattice with 
weak disorder. The atoms are assumed to be initially prepared in a superfluid state 
with non-zero quasimomentum and, hence, may carry matter current. It is found that 
the atomic current persists in time for a low value of the quasimomentum but decays 
exponentially for a high (around one quater of the Brillouin zone) quasimomentum. 
The explanation is given in terms of low- and high-energy spectra of the Bose- 
Hubbard model, which we describe using the Bogoliubov and random matrix theories, 
respectively. 
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1. Introduction 

Ultracold atoms in optical lattices constitute an intense research activity both in 
experimental and theoretical physics. Up to now this system has mostly been used for 
modelling the fundamental Hamiltonians of solid state theory (see, PQE], for example) 
where the number of particles is macroscopically large. However, the recent progress 
with manipulating a countable number of atoms Hj makes it possible to build a 
system of arbitrary size, ranging from microscopic to macroscopic. In this border region 
between microscopic and macroscopic one has to deal with a finite number of atoms 
which, on the one hand, is too large to use the single-particle approach but, on the 
other hand, is too small to justify the thermodynamic limit. In the present work we 
theoretically analyse one of these problem related to superfluidity of a few (N ~ 10) 
Bose atoms in a ring optical lattice |H] with a few (L ~ 10) sites. 

It should be stressed in the very beginning that, currently, there are two different 
definitions of superfluidity in the physics literature. One definition is based on the 
system's responce to a phase twist. With respect to Bose atoms in a lattice this approach 
is discussed, in particular, in Ref. [Hj, and a method of how one can realize the twisted 
boundary conditions in a laboratory experiment is suggested in Ref. [Sj. The other 
definition originates in the Landau criterion of superfluidity and involves a responce 
of a superfluid flow to 'wall roughness' [HE]- In this work we try to reconcile both 
approaches. Specifically, we address the following problem. Assume that we have N 
Bose atoms in a ring lattice with L sites in a superfluid state with given quasimomentum 
k = 2nk/L: 



We are interested in the time evolution of this state (which we also shall refer to as 
the supercurrent state) in the presence of a weak scattering potential and atom-atom 
interactions. 

We note that for a BEC of atoms (N ^> 1) the problem of superfluid atomic current 
has been considered in a large number of papers (see Refs. PI ITUl ITI] IT21 IT31 IT^l ITo] ITo] . 
to cite few of them). The starting point of all these studies is the mean-field approach, 
which is sometimes rectified by taking into account the quantum fluctuations [T5| I16j. 
The mean-field theory predicts a destruction of the supercurrent as soon as the 
quasimomentum exceeds one quater of the reciprocal lattice constant (n > n/2 in the 
notations used). In order to justify the mean-field approach in a ID lattice the mean 
number of atoms per one site should be much larger than unity. As stated above, in the 
present work we focus on the opposite limit N/L ~ 1, where the mean- field approach is 
not applicable. For this reason we treat cold atoms in an optical lattice from a different 
viewpoint, in a sense closer to quantum optics than to condensed matter physics. 

The paper essentially consists of two parts, - in the first part (Sec. EJ), after a 
brief preliminary analysis, we report the results of numerical simulations of the system 
dynamics, and in the second part (Sec. El and Sec. EJ) we explain the observed regimes 




(1) 
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in terms of the energy spectrum of the system. The main results are summarized in the 
concluding Sec. El 

2. Supercurrent dynamics 

Before proceeding with numerical simulations, we shall briefly discuss possible regimes 
for the atomic current. 

2.1. Preliminary analysis 

Let us first consider the single-particle problem. In the tight-binding approximation the 
Hamiltonian of the system reads 

ff = ~£(|J + i><i| + *-c.) + £W<J|. (2) 
z i i 

where |/) are the Wannier functions, J the hopping matrix element, and V\ the random 
scattering potential. In what follows, to be concrete, we shall consider < V/ < e 
with e C J. The operator V = J2i ^ 10(^1 couples the degenerate states with opposite 
quasimomentum, resulting in new eigenstates \k cs ) = (\k) ± | — k))/\/2 with an energy 
splitting \E C -E S \ = 2\(k\V\ - k)\ = 2|V(2fc)|, where 

K(*) = iEvi«p(^«)~f 

is the Fourier transform of VJ. Thus, in course of time, an atom in a ring will periodically 
change its momentum to the opposite one with the frequency Q e ~ e/hL. It is worth 
of stressing that this periodic dynamics is exclusively due to the fmiteness of L and the 
assumed condition e < J, which means that V couples only the degenerate states of 
the unperturbed system, f 

Next we consider the multi-particle case, 

h = ~ £ (sf +1 a* + h.c) +^X>(ri,-i)+x; vm , (3) 
z i z i i 

where a\ and di are the bosonic creation and annihilation operator, hi = ajdi, 
and U is the on-site interaction energy. Using the canonical transformation, bk = 
(l/VL) J2i exp(i2irkl / L)di, it is convenient to present the Hamiltonian © in the form 




(4) 

where S(k) = 1 if A: is a multiple of L, and 5(k) = otherwise. For U = and 
e = the multi-particle eigenstates of the system are the quasimomentum Fock 

| In terms of Anderson's localization theory the above conditions mean that the Anderson localization 
length is much larger than the system size. 
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states |n) = \n , n\, . . . , Ul-i), where Y,k n k — N. Our state of interest corresponds 
to \k) = | . . . , 0, Nk, 0, . . .), where all atoms have one and the same quasimomentum. 
Similar to the single-particle case, the random potential couples this state to the 
supercurrent state with the opposite quasimomentum | — k) = \ . . . , 0, N^, 0, . . .), 
k' = modi (—&;). However, now the coupling is indirect and involves the intermidiate 
states \n{m)) = \ . . . , (N — m) k , . . . , m^, . . .), as it immediately follows from the explicit 
form of the scattering potential in the momentum representation. Thus the time 
evolution of the state \k) is defined by the following (N + 1) x (iV + 1) matrix, 



+ y/(N - m)(m + 1) [V(2k)5 m+ i, m > + V*(2k)S 

m,m'+l ] ,(5) 

where E m — E K = — JN cos k are the degenerate energies of the states |/t(m)) and the 
next terms the transition matrix elements {K(m)\V\K(m')) . The spectrum of the matrix 
(0) is equidistant with the level spacing 2|V(2/c)|. Thus we have reproduced the result 
of the single-particle analysis, where the time evolution of the system is periodic with 
the frequency Q e ~ e/TiL. 

Now we switch on the interaction. Then the intermidiate states \n{m)) aquire 
energy shifts E m = E m (U), which appear to be m-dependent. Using the first order 
perturbation theory we obtain 

E m = E K + ^(^m^^blhlh^bkjih + k 2 - k 3 - fc 4 )|/t(m)) 

= E K + — [(N - m) (N - m - 1) + m{m - 1) + 4(iV - m)m] 
UN 2 U 

~ E K + — — + -m(N - m) . (6) 

ZLj Lj 

Due to the mismatch of the energy levels E m , the supercurrent states \k) = | /t(0) ) and 
| — k) = |^(^)) become effectively decoupled and, hence, the supercurrent should persist 
in time. 

At this point we would like to note the analogy of the problem discussed with that 
for a BEC in double well potential . Drawing this analogy further we can estimate 
the minimal U m i n required for stabilization of the suppercurrent as 

U mm « 8e/N . (7) 

For the sake of completeness we present a derivation of the estimate ((7|) in the next 
subsection (which can be safely skipped if a reader is familiar with the subject). 



2.2. Semiclassical approach 

The standart method of treating the system (jHHHI) consist of mapping it onto an effective 
classical system (terms proportional to the identity matrix are omitted), 

H eff = gl(l - I) + 2\V\y/l(l - I) cos9 , g = UN/L, (8) 

followed by a semiclassical quantization, where 1/N plays the role of Planck's constant. 
The phase portrait of the system (jHJ) is shown in Fig.^for g/|U| = 1 and = 10. It 




Figure 1. Phase portrait of the effective system JHJ) for g/\V\ = 1 (left panel) and 
g/\V\ = 10 (right panel). 



is seen that when g is increased the phase portrait becomes similar to that of classical 
pendulum, with a separatrix separating the librational and rotational regimes. The 
maximal and minimal values of the classical action / along the separatrix are given by 
I* ps 1/2 ± J\V\/2g. The quantum states, associated with /, are decoupled only if they 
lie above the separatrix. Then, by requiring \I* — 1/2 1 < 1/2 and noting that \V\ ~ e/L, 
we come to the estimate (J7|). 

Needless to say, the semiclassical approach described above requires 1/N <C 1 and 
is not accurate for small N. Nevetheless, even for N ~ 10 the spectrum of the matrix 
A can be well understood in terms of the effective system (jHJ). For the purpose of future 
reference, the right panel in Fig. El shows the numerical solution of the matrix eigenvalue 
problem for N = 7, L = 9, and \ V\ = 0.0168. In particular, at U — 0.2 J one can identify 
the first four top levels with the phase trajectories below the separatrix, next two levels 
with trajectories around the separatrix, and the last two almost degenerate levels with 
trajectories well above the separatrix. 

2.3. Numerical results 

The above conclusion about the persistent current relies on the applicability of a 
perturbative approach. Formally this means that the supercurrent state \k), as well 
as the intermidiate states |/c(m)), have to be approximate eigenstates of the system 
at e = 0. This imposes the upper boundary U max on the interaction constant, which 
appears to depend on the quasimomentum k. Indeed, the state \k) with the energy 
E K — J N cosn is an approximate eigenstate of the system only if U is smaller than 
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Figure 2. Spectrum of the matrix (l- r >l(i[l as function of the interaction constant U. 
(The energy is measured in units of J, Eq — — J N + U N 2 /2L .) Parameters are N = 7, 
L = 9, k = 2tt/L, and |V(2)| = (left panel) and \V(2)\ = 0.0168 (right panel). 



the characteristic energy gap separating it from the other energy states, coupled to \k) 
by interaction. As the first guess one can set this gap to the mean level spacing, given 
by the inverse density of state AE = l/f(E). It is easy to show that for Uf J < 1 the 
density of states of (jHJ) is given by the Gaussian distribution (see Fig. HI below) 



f(E) 



M 



2na 



exp 



{E-Ef 
2o^~ 



(9) 



where M = (N + M - 1)!/JV!(M - 1)! is the dimension of the Hilbert space, a ~ JVN 
and E ~ UN 2 /L + eN/2. Thus the characteristic gap for a supercurrent state, which 
belongs to the central part of the spectrum (i.e., for k ~ vr/2), is essentially smaller than 
that for a supercurrent state with low quasimomentum k 7r/2. As a consequence, 
U m ax f° r the supercurrent state with k ~ it/2 may be smaller than U m i n . In the other 
words, the perturbative approach of Sec. 2.1 (where we used first order perturbation 
theory to find corrections to the eigenenergies of states |ft(m))) breaks down before the 
stabilization of the supercurrent is achieved. 

Figure El compares the dynamics of iV = 7 atoms in a lattice with L = 9 sites, 
which were initially prepared in the supercurrent state (0) with high, k = 6n/L, and 
low, K = 2ir/L, quasimomentum. The normalized mean momentum of the atoms, 
pit) = A r_1 Im[(\I / (t)| J2i a] + iai\ty(t))], is depicted. It is seen in the upper panel of Fig. El 
that in the former case of high quasimomentum the oscillatory behaviour of p(t) changes 
to irreversible decay as the interaction constant is varied from U = 0.02 to U — 0.2. 
(From now on we set J = 1, i.e., energy is measured in the units of J and time in the 
units T = 2irh/J.) Further increase of the interaction constant (results are not shown) 




1 10 20 30 40 

t/T 

Figure 3. The mean momentum of N = 7 atoms in a lattice with L = 9 sites. The 
magnitude of the scattering potential e = 0.2. The interaction constant U = 0.02 
(dash-dotted lines), U = 0.1 (dashed lines), and U — 0.2 (solid lines). The atoms 
are initially prepared in the supercurrent state Q with n — 6n/L (upper panel) and 
k = lit I L (lower panel). 



is reflected in even faster decay of p(t). This should be contrasted with the case of low 
quasimomentum (lower panel), where the current oscillations at U — 0.02 change to 
persistent current at U — 0.2. Here further increase of U leaves the system dynamics 
qualitatively unchanged at least till U = 1. The displayed numerical results suggest 
that the perturbative approach of Sec. 2.1 works for k = 2ix / L but does not work for 
k > 2tt/L. We shall come back to this point later on in Sec. 4.3. 

It is worth stressing that through the paper we consider a single real- 
ization for the random potential (i.e., no average over disorder). Specifically 
to the considered lattice of L = 9 sites, the random entries are V\ = 
e(0.80, 0.59, 0.06, 0.18, 0.97, 0.31, 0.67, 0.78, 0.49). The Fourier transform of this sequence 
gives \V(2k) \ = 0.084e and |Vi2Jfe)| = 0.144e for jfe = 1 and k = 3, respectively. 

3. High-energy spectrum 

To get a better insight in the physics of the discussed phenomena we shall discuss the 
displayed in Fig. El results in terms of the energy spectrum of the system (jSJ). We begin 
with the case of a high quasimomentum which, as mentioned above, refers to the central 
part of the spectrum. 
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3.1. Spectral statistics 

We have found that in the case of high initial quasimomentum a transition from 
oscillatory dynamics to irreversible decay is associated with the transition to chaos 
in the Bose-Hubbard model. Following Ref. ^H], we shall monitor this transition by 
analysing the distribution of distances between the neighbouring levels, normalised to 
the mean level spacing: s = (E n+ i —E n )/AE = (E n+ % — E n )f[(E n+ i + E n )/2]. It should 
be stressed that the presence of random potential in the Hamiltonian @ alone does not 
yet induced chaos in the system. The only consequence of a weak disorder (relevant to 
the spectral statistics) is that it breaks the translational symmetry and, hence, we need 
not worry about decomposition of the energy spectrum into the independent subsets 
(labeled, in the absence of a random potential, by total quasimomentum of the atoms 



The results of the statistical analysis of the high-energy spectrum are presented 
in Fig. |U The dash-dotted and dashed lines in panel (c) correspond to the integrated 
distribution, I(s) = J S P(s')ds', for the Poisson statistics, 



which is typical for a generic integrable system, and the Wigner-Dyson statistics, 



typical for non-integrable systems. These distributions reflect the different character 
of the parametric dependence of the energy levels E n = E n (X) on some parameter in 
the Hamiltonian (A = U in our case). Namely, in the integrable case the energy levels 
may cross and, hence, one finds an arbitrary small s. On the contrary, if the system is 
non-integrable, the energy levels show avoided crossings and probability of finding small 
s tends to zero. 

The panel (a) in Fig. 0] shows the density of states f(E) for U = 0.02, where 
only the data from the central part of the spectrum (marked by the inverse parabola) 
were used for the statistical analysis. It is seen in the lower panel that for U = 0.02 
the level spacing distribution follows the Poisson statistics. Thus for this value of the 
interaction constant the system should be classified as integrable, which is consistent 
with the periodic dynamics of the mean momentum in Fig. Of a). The panel (b) in 
Fig. |3] shows the density of states for U = 0.2. Apart from an uniform shift of the 
spectrum to positive values, no qualitative change in f(E) is observed. However, we do 
observe a qualitative change in the level spacing distribution. Now it reliably follows the 
Wigner-Dyson statistics which, as mentioned above, is a hallmark of quantum chaos. 

3.2. Local density of states 

The spectral statistics is only one (and, in fact, rather poor) characteristic of the system. 
In particular, the level spacing distribution remains unchanged (Wigner-Dyson) in the 
interval 0.2 < U < 1, although the decay rate of the supercurrent changes with U. One 



[19). 



P{s) = exp(-s) , 



(10) 




(11) 





Figure 4. (a,b) - density of states of N = 7 atoms in a lattice with L — 9 sites for 
the interaction constant U = 0.02 and U = 0.2, respectively. The magnitude of the 
scattering potential e = 0.2. (c) - integrated level spacing distributions for the central 
part of the spectrum. 



gets more information about the system by studying its eigenf unctions. To this end we 
introduce a quantity R(m,n), 

R(m,n) = \{* m (U')\V n (U))\ 2 , (12) 

closely related to the so-called local density of states. § In Eq. (JT2"j) \ty n (U)) are the 
eigenfunctions of the Hamiltonian © calculated for a given U and ordered according 
to their energies. In what follows we shall fix U' = 0.02 while U will be scanned in the 
interval 0.2 < U < 1. Since for U = 0.02 the system is integrable, the matrix (|12|) can 
be alternatively viewed as the matrix of the expansion coefficients of the chaotic states 
\^ n (U)) over 'regular basis' \m) = \$f m (U = 0.02)). 

The characteristic structure of the matrix (|12jl is shown in Fig. El for U = 0.2. It 
is seen that that R is a banded matrix with strongly fluctuating matrix elements. The 
mean values of the elements across the main diagonal, 

j M/2 

R(Am) = — R{m,m + Am) , ^ii(Am) = l, (13) 

m=-M/2 Am 

are shown in Fig. Elon linear and logarithmic scales. (Here, as in the spectrum analysis, 
we consider an energy window of the order of unity in the central part of the spectrum.) 

§ The local density of states is defined as R(m, E) = J2 n R( m i n)5(E — E n ). 
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It is seen that R(Am) converges to the Lorentzian, || 

R(Am) = -r T J 2n ^ l . (14) 
v ; (Am) 2 + r 2 /4 v ; 

We note, in passing, that a similar result is reported in the recent paper [2T| devoted to 
the spectral properties of the three-site Bose-Hubbard model. 

The distribution (|14j) . also known as the Breit-Wigner formula, implies the 
exponential decay of the supercurrent state. Indeed, considering the overlap integral 
(k | «;(£)), one has 

(ft|exp (--Htj | «) = (K\m)(m\^ n )exp (~E n t\ (^ n \m'){m'\K) 

m,m'n 



where we substitute the exact energy levels E n by their approximate positions, E n 



E K +AEAm. (This approximation obviously holds till time t% ~ fi/AE = hf(Ek), which 
increases exponentially with the system size.) Substituting here R(Am) from Eq. (fTlj) 
we have (n\K(t)) = exp(—TAEt/h). We found that the width T grows approximately 
quadratically with U in the interval 0.2 < U < 1. 

|| The distribution (|14l) is typical for the banded random matrices |2U) . It is interesting to note in this 
connection that for the N/ L ~ 1 neither matrix of the Hamiltonian (J3J) nor that of the Hamiltonian 
ijljl are banded. It is an open problem in the random matrix theory to extend the results of |20| to the 
present case of very sparse but not banded matrices. 
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Figure 6. Mean values of the matrix elements across the main diagonal in the central 
part of the matrix. The solid line is the best fit by the Breit-Wigner formula l|14|l . 

4. Low-energy spectrum 

We turn to the case of low quasimomentum. For small k the energy of the supercurrent 
state falls into the low-energy tail of the density of states (JHJ), where the random 
matrix approach is not applicable. On the other hand, the low-energy spectrum of the 
interacting Bose atoms is believed to be described by the Bogoliubov theory. For this 
reason we review the Bogoliubov approach for a finite size system. Through the section, 
if not stated otherwise, we assume the homogeneous case e = 0. 

4-1. Bogoliubov approach 

As an intermidiate step, let us show that the Bogoliubov approach amounts to the 
following two assumptions, (i) The low energy eigenstates of interacting Bose atoms are 
given by a linear superposition of the quasimomentum Fock states, where n atoms 
have quasimomentum k, n quasimomentum — k, and the rest N — 2n have zero 
quasimomentum [IB], i.e., 

N/2 

= H c n\N - 2n, ... ,n k , , 0) , k' = mod L (-k) . (15) 

n=0 

(ii) The number of atoms with k ^ is small compared to the number of atoms with 
zero quasimomentum, i.e., only the coefficients c n with n <C N/2 have non-negligible 
values. (This condition is automatically satisfied if one assumes the thermodynamic 
limit iV -»• oo, U -> 0, UN = const.) 
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The analysis goes as follows. Substituting the wave function (JT5J) in the eigenvalue 
equation with the Hamiltonian (£Q), we get a system of linear equations for the coefficients 

(2E k n + ^-a n )c n + ^(6 n c„_i + b n+1 c n+1 ) = Ec n , (16) 

where 

a n = 2nN - 3n 2 + n + N(N -l)/2« n{2N - 3n) + N 2 /2 , 

b n = (n + l)y/(N -2n)(N - 2n - 1) « n(JV - 2n) , 

and = J[l — cos(27r fc/L)] is the single-particle excitation energy (should not be 
missmatched with the energy of the supercurrent state, E K = — JNcosn). Assuming 
the thermodynamic limit, Eq. (|16j) simplifies to 

2(E k + g)nc n + gnc n ^ + g(n + l)c n+1 = Ec n , g = NU/L. (17) 

Next, introducing the generating function, 

1 



V Z7T 



And 



n=—oo 



we present the system of linear equations ()17j) as a differential equation on the function 

S [ne* e + 2(1 + e)n + e~ ie n] <$>(6) = E$(0) , (18) 
where n = —id/ 89 and e = E k /g. The general solution of ()18j) reads 

$(0) = Cexp f i f E/9 n - ^ d-&\ . (19) 
v ; y \ Jo 2cos$ + 2 + e J K ' 

Finally, requiring Q(8 + 27r) = $(0) and calculating the relevant integral, 

1 dd 1 



2n Jo 2cos6» + 2 + £: 2^(1 + e) 2 - 1 
we get the equidistant spectrum with the transition frequency 



u k = 2pgE k + E 2 k . (20) 

The result (|20|) reproduces the famous Bogoliubov equation for the quasiparticle 
excitations of the Bogoliubov vacuum. 

4-2. Bogoliubov spectrum 

In the previous subsection we have considered an excitation of the given quasimomentum 
state, with the single-particle excitation energy E k = J(l — cos k). To include the other 
quasimomentum states, the ansatz (fTH|) should be generalized to 

|tf> =J2c n \N-2^2n k ,n 1 ,n 2 ,...) , (21) 

n k 

% To be regorous, Eq. Ijl8|l is not strictly equivalent to Eq. (|1T|) in the sense that it also has solutions 
with negative E. 
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where n = (n\, . . . , n L / 2 ). Substituting (|2Tj) in the stationary Schrodinger equation with 
the Hamiltonian (jljl. we obtain a system of rather complex equations on the coefficients 
c n , which can be solved analytically only in the thermodynamic limit. In this limit, as 
it is easy to show, the whole eigenvalue problem factorizes to L/2 eigenvalue problems 
of the form (j!7)) and, hence, the whole spectrum is given by the direct sum of L/2 linear 
spectra. 

A remark about the total quasimomentum, which is a global symmetry of the 
system in the absence of random potential, is in turn. The substitution (fTH|) corresponds 
to zero total quasimomentum. To get non-zero values of the total quasimomentum, one 
should use a slightly different ansatz, 

N/2 

= J2 °n\N - 2n - m, . . . , n k + m, . . . , n' k , . . . , 0) , (22) 

n=0 

k! = modii—k) , m — 1, . . . , L — 1 . 

Ansatz (f2*2*j) leads to the eigenvalue equation of the form (fTKj) but with different 
coefficients a n and b n . In particular, considering the thermodynamic limit, Eq. (fTTj) 

changes to 

2(E k + g)(n + m)c n + g\Jn{n + m)c n _i + g\J (n + l)(n + m + l)c n+ i = Ec n . (23) 

We note, in passing, that if the spectra associated with different single-particle excitation 
energy and different total quasimomentum are superimposed, 

E = Y,{E {Km \g)} , (24) 

k,m 

one typically finds a multiple degeneracy of the levels at U = (see Fig. El below). 

It is interesting to compare the discussed Bogoliubov spectrum of an infinite system 
with the low-energy spectrum of a finite system. For this reason we calculate numerically 
the spectrum of iV = 25 atoms in a lattice with L = 5 sites. + A few first levels of 
this system are depicted in the left panel of Fig. [7| where, to facilitate a comparison, 
we substract the energy Eq = —JN + UN(N — 1)/2L and rescale energy axis on the 
basis of the frequency uj\ = Wi(U). The Bogoliubov spectrum, calculated by using 
Eqs. fl23H24|) . is depicted in the right panel of Fig. Both similarities and differences 
are evident. The first two levels are seen to coincide in the whole interval < U < 1. 
On the other hand, multiple degeneracy of the levels around f7 = 0.13 in the right panel 
is removed in the left panel. This is, in fact, not surprising. Indeed, let us consider 
the lowest group of levels, showing the degeneracy. These levels are associated with the 
qusimomentum Fock states |0, 2, 21, 2, 0), |0, 2, 22, 0, 1), |1, 0, 22, 2, 0), and |1, 0, 23, 0, 1). 
(Here we use a different notation for the Fock states, corresponding to the Brillouin zone 
—Ti < k < n.) These states belong to different spectra, labeled by k and m in Eq. (|2HJ), 
and are decoupled within the Bogolubov approach. However, for the considered finite 
system these states are coupled by interaction, where the coupling matrix elements are 
of the order of g/N. 

+ For L = 5 there are two different frequencies Wk- In this sense, L = 5 is the simplest generic case to 
discuss the Bogoliubov spectrum. 
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Figure 7. Left panel: first few energy levels of N = 25 atoms in a lattice with L = 5 
sites, as function of the on-site interaction constant. (Only the levels corresponding 
to zero total quasimomentum are shown.) Right panel: first few energy levels of the 
Bogoliubov spectrum, calculated on the basis of Eq. and Eq. (PUl . 



4-3. Persistent current 

In this subsection we critically review the result of Sec. El about the persistent current, 
carefully checking validity of the perturbative approach. Left panel in Fig. |H] shows the 
low-energy levels for N = 7 and L = 9 (whole spectrum is shown, i.e., no symmetry 
selection according to the total quasimomentum). Remarkably, even for such a small 
number of atoms one still has a qualitative agreement with the Bogoliubov spectrum. 
Our states of interest in Fig. IHlare the supercurrent and the intermidiate states |«;(m)), 
which originate from the point marked by an asterisk. It is seen, by comparing with 
Fig. Efa), that (i) the splitting between these levels well matches Eq. © and (ii) the 
coupling of these states to the other states of the system is negligible, which is indicated 
by absence of the avoided crossings. In addition to the case e = 0, right panel in Fig. El 
shows the spectrum of the atoms in the presence of a weak scattering potential [should 
be compared with Fig. Efb)]. Again, no avoided crossings with the other levels are seen. 
Hence, the approach of Sec. 2 is well justified. 

The above visual analysis of the spectrum can be made quantitative by considering 
the overlap of the supercurrent state \k) with the exact eigenstates, 

Q(U) = m a x(\(^ n (U))\ 2 ) . (25) 

71=1 V / 

For k = the quantity (J25J) is obviously maximazed by the ground state \^o(U)), which 
is expected to coincide with the Bogoliubov state ground. The solid line in Fig. 
shows the overlap of the state \k = 0) with the ground Bogoliubov state. A monotonic 
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Figure 8. Low-energy levels of N = 7 atoms in a lattice with L — 9 sites for magnitude 
of the scattering potential e = (left) and e = 0.2 (right). The state with supercurrent 
(k = 2tt/L) corresponds to the lowest level in the group of levels marked by an asterisk. 
Open circles corresponds to the quasiparticlc energy h,u>\. 

decrease of Q = Q(U), seen in the figure, is due to population of the single-particle 
quasimomentum states with k^0, and is often referred to as the Bogoliubov depletion 
of the BEC. Additionally, the dashed and dash-dotted lines in Fig. |Hl depict the overlap 
of the states \k = 2n/L) and \k = Att/L) with the Bogoliubov states, originating 
from these supercurrent states, which we calculate by substituting the single-particle 
excitation energy E k = J[l — cos(27ik/L)} in Eq. (JT7J) by 

E k = 0.5 J[cos(k + 2nk/L) + cos(k - 2nk/L) - 2 cos k) ~ (2nk/L) 2 cos k . (26) 

Finally, the series of dots correspond to the quantity (|25jl . It is seen that for k = 
dots perfectly follow the solid continuous line. Thus the ground state of the system 
is indeed well approximated by the Bogoliubov state. With exception of two narrow 
avoided crossings this is also the case for the state of our interest k = 2ir/L. However, for 
higher initial quasimomentum k = 47r/L, the Bogoliubov state is seen to be completely 
destroyed by the large number of avoided crossings. 

5. Conclusions 

Within the formalism of the Bose-Hubbard model we have considered time evolution 
of the atomic supercurrent in a ring optical lattice with weak on-site disorder. For 
vanishing atom-atom interactions, weak disorder induces Rabi oscillations of the atomic 
current, where the atoms periodically change their velocity to the opposite one. For 
non-vanishing atom-atom interactions, the supercurrent dynamics depend crucially on 
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Figure 9. Overlap of the states \k = 0) (solid line), |re — 2tt/L) (dashed line), 
and | re = Att/L) (dash-dotted line) with the Bogoliubov states, originating from these 
supercurrent states. Dots (guided by the solid lines) show the quantity l)25[l. calculated 
for these three values of the initial quasimomentum. The system parameters are the 
same as in Fig. |Hfa). 

the initial quasimomentum k (i.e., the initial velocity of the atoms). Namely, for a 
high quasimomentum k ~ ir/2 the supercurrent exponentially decays as the interaction 
constant U exceeds some critical value, while for a low quasimomentum k <C 7r/2 the 
oscillatory behaviour of the suppercurrent changes to a persistent current. 

The explanation for these effects is found in the structure of low- and high-energy 
spectra of the Bose-Hubbard model. It is shown that the low-energy spectrum of 
the system is regular, and the positions of the energy levels can be found by using 
a Bogoliubov approach. In contrast, the high-energy spectrum shows a transition from 
regular to a chaotic one if U exceeds its critical value. Using the results of the random 
matrix theory, we show that this transition is reflected in the exponential decay of the 
supercurrent with the decay constant proportional to U 2 . 

This work was supported by Deutsche Forshungsgemeinschaft within the SPP1116 
program. 
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